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ABSTRACT 


Measurement  of  the  absolute  amplitudes  of  a  seismic  signal  requires  accurate  knowledge  of  the 
seismometer-digitizer  response.  When  the  seismometer  and  digitizer  characteristics  are  known,  deconvolution  of  the 
digital  system  response  is  relatively  straightforward.  However,  even  with  a  known  response,  problems  in  timing  and 
incorrect  polarities  can  often  contaminate  any  calibration  parameters  derived  from  these  data.  When  either  one  or 
both  components  of  the  response  for  a  seismometer  system  are  unknown,  the  data  are  often  not  considered  suitable 
for  seismological  tasks  that  require  absolute  amplitude  measurements. 

We  have  developed  a  set  of  software  tools  that  can  recover  the  sensor/digitizer  response  function  from  raw  data 
when  the  response  is  either  unknown  or  incorrectly  known.  These  software  tools  include  SACPSD,  which  estimates 
seismic  noise  power  spectral  densities,  and  NOISETRAN,  which  generates  a  pseudo-amplitude  response  (PAR)  for 
a  seismic  station,  based  on  background  noise  at  a  reference  station. 

During  the  first  year  of  our  Phase  II  project,  we  have  focused  on  improving  NOISETRAN  and  validating  the  method 
on  multiple  datasets.  First,  we  estimate  the  PAR  for  a  short-period  seismometer  (L-4C  3D)  using  a  broadband  station 
located  40  km  away.  The  PAR  was  within  a  factor  of  2  of  the  true  response  for  the  seismometer  from  periods 
between  0.1  and  3  Hz.  Next  we  applied  the  method  to  USArray  data  from  2009-2011.  Application  to  the  400+ 
stations  of  the  Transportable  Array  (TA)  deployed  in  the  central  United  States  and  Rocky  Mountains  allowed  us  to 
understand  spatial  variations  in  the  response  error,  which  is  the  difference  between  the  actual  response  and  the  PAR. 
For  stations  within  200  km  of  the  reference  station,  the  response  error  for  all  frequencies  (0.01  to  10  Hz)  is  typically 
within  less  than  ±0.3  log  units  (e.g.,  factor  of  ±2  difference  between  the  actual  and  pseudo-response).  The  scatter  in 
the  error  estimates  for  all  frequencies  increases  to  ±1.0  log  units  (e.g.,  order  of  magnitude)  at  great  distances; 
however,  much  of  that  error  is  at  frequencies  above  1  Hz,  where  local  cultural  noise  source  differences  between  the 
reference  and  candidate  stations  are  magnified.  We  also  see  geologic  and  population  effects  distributed  in  the 
response  error  maps.  Lastly,  we  used  the  PAR  files  estimated  from  NOISETRAN  to  correct  the  TA  data  from  the 
1 1  March  2011  Tohuku  earthquake.  The  error  between  maximum  P-wave  amplitudes  estimated  using  the  PAR  files 
and  the  actual  responses  ranges  from  -0.3  log  units  (for  frequencies  below  0.5  Hz)  to  -1.0  log  units  (for  frequencies 
above  0.5  Hz).  The  response  error  is  consistently  negative,  meaning  that  the  PAR  files  are  overcorrecting  the 
P-wave  amplitudes. 

During  the  final  year  of  the  project,  we  will  package  NOISETRAN  and  additional  seismic  tools  into  an  easy-to-use 
graphical  interface,  conduct  alpha  testing  of  the  toolbox,  and  then  conduct  beta  testing  at  Los  Alamos  National 
Laboratory. 
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OBJECTIVES 


We  have  developed  a  seismie  software  toolbox  with  the  objeetive  of  estimating,  assessing,  and  ehronieling  the 
seismometer  and  digitizer  response  of  a  seismie  station  when  the  aetual  response  is  unknown.  During  the  first  year 
of  this  Phase  II  projeet,  we  initiated  development  of  an  easy-to-use  graphieal  user  interfaee  to  assist  the  user  in 
applying  these  tools  to  iteratively  reeover  the  response  funetion  of  an  unknown  instrument.  We  have  foeused  on 
eontinued  development  and  validation  of  NOISETRAN,  whieh  estimates  the  seismie  response  of  a  station  by 
eomparing  its  baekground  noise  and  reeorded  signals  to  noise  and  signals  from  nearby  stations  with  known 
responses.  In  this  paper,  we  provide  an  example  of  the  NOISETRAN  method  using  seismie  data  from  Israel, 
estimate  error  and  uneertainty  in  the  method  using  the  USArray  stations,  and  then  apply  the  NOISETRAN-estimated 
responses  to  reeordings  of  the  2011  Tohoku  earthquake. 

RESEARCH  ACCOMPl  ISHED 

NOISETRAN:  Noise  Transfer  Functions  for  Pseudo-Response  Estimation 

Important  ealibration  datasets  may  be  obtained  in  whieh  no  metadata  (e.g.,  station  information,  sueh  as  loeation, 
instrument  type,  response,  and  gains)  aeeompany  the  waveforms.  We  are  examining  methods  that  attempt  to 
determine  a  rough  order-of-magnitude  (ROM)  estimate  of  the  seismometer-digitizer  response  by  eonsidering  the 
baekground  seismie  noise  at  a  nearby  seismie  station  with  a  known  response  (e.g.,  the  “referenee”  station).  We  refer 
to  this  method  as  NOISETRAN,  as  we  attempt  to  estimate  a  transfer  function  T{k>)  to  eonvert  the  baekground 
seismie  noise  at  a  eandidate  station  (R X®))  with  unknown  response  (heneeforth  referred  to  as  the  “eandidate” 
station)  to  the  baekground  noise  levels  at  a  referenee  station  R2(ro).  The  transfer  funetion  ean  be  estimated  using: 

Ri{a))^R2{(£))*T{(i>)+er((i>).  (1) 

The  response  error  (e^)  in  this  proeedure  is  due  mainly  to  differenees  in  the  noise  fields  between  the  two  sites,  whieh 
is  frequeney  (®)  dependent.  Onee  the  ROM  system  response  Rfto)  has  been  reeovered,  we  ean  then  use  waveform 
modeling  teehniques  to  further  refine  the  system  response. 

Our  assumption  for  attempting  the  NOISETRAN  teehnique  is  that  the  aetual  noise  field  at  eertain  periods  should  be 
the  same  or  similar  for  a  network  eovering  a  relatively  small  geographieal  area.  MeNamara  and  Buland  (2004) 
showed  the  geographieal  distribution  of  noise  in  three  different  period  bands  (e.g.,  0.125-0.0625  see,  8—4  see,  and 
64-32  see)  in  North  Ameriea.  They  observed  the  strongest  geographieal  variations  at  periods  smaller  than  1  see  due 
to  eultural  noise,  espeeially  near  large  population  eenters.  Proximity  to  oeeans  dominates  the  variations  in  the 
mieroseism  band  (8-4  see).  They  note  that  at  longer  periods  (64-32  see)  the  geographie  variations  deerease  and  that 
the  responses  are  largely  a  funetion  of  the  installation  eharaeteristies  (e.g.,  borehole  vs.  surfaee  installation)  instead 
of  the  aetual  seismie  noise  eharaeteristies.  Their  results  suggest  that  our  best  results  for  NOISETRAN  will  most 
likely  oeeur  at  periods  longer  than  1  see. 

Applying  NOISETRAN  to  a  Seismometer  in  Israel 

During  Phase  I  of  this  projeet,  we  applied  NOISETRAN  only  to  broadband  seismie  stations,  whieh  have  typieal 
responses  that  are  flat  aeross  a  wide  range  of  frequeneies  (e.g.,  0.01  to  50  Hz).  The  first  test  for  the  Phase  II  projeet 
was  to  determine  the  applieability  of  the  NOISETRAN  method  on  instruments  typieally  referred  to  as  “short-period” 
seismometers  (e.g.,  L4,  L22,  S6000)  or  “high-frequeney”  geophones  (L28,  ete).  These  instruments  are  typieally  flat 
to  veloeity  at  frequeneies  above  1,  2,  or  4.5  Hz,  but  have  a  steep  roll-off  at  lower  frequeneies. 

In  January  2011,  Weston  Geophysieal  Corporation  (WGC)  partieipated  in  the  Sayarim,  Israel,  Ground  Truth  Zero 
(GTO)  Calibration  explosions  (Gitterman  and  Hofstetter,  2011),  whieh  ineluded  a  102 -ton  surfaee  explosion  used  for 
ealibrating  the  International  Monitoring  System  (IMS)  network.  WGC  deployed  six  (6)  Sereel  L-4C  3D  veloeity 
seismometers  within  10  km  of  the  explosions.  Although  we  already  knew  the  response  for  the  L-4C  3D 
instrument — ^the  seismometer  is  fiat  to  veloeity  above  1  Hz — we  used  NOISETRAN  to  reeover  the  response  of  one 
of  the  stations  (NS3)  deployed  for  the  Sayarim  GTO  experiments.  We  seleeted  a  sample  of  nighttime  noise  for  the 
eandidate  station  for  input  into  NOISETRAN.  The  next  step  in  NOISETRAN  is  to  seleet  a  referenee  station.  The 
eonsiderations  for  ehoosing  a  referenee  station  should  inelude  (1)  available  data  to  establish  a  baekground  noise 
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model  for  the  reference  stations,  (2)  proximity  to  the  stations  of  interest,  and  (3)  whenever  possible,  similar  expected 
background  noise  levels  between  the  reference  and  unknown  stations.  The  station  that  best  matched  these 
requirements  for  our  candidate  station  was  EIL  (Eilat,  Israel).  We  downloaded  data  from  EIL  from  the  Incorporated 
Research  Institutions  in  Seismology  (IRIS)  for  the  same  time  period  as  our  candidate  station  noise  sample. 

Figure  1  compares  the  power  spectral  densities  (PSDs)  for  the  candidate  and  reference  station  after  the  response  fde 
was  repaired.  We  estimate  the  transfer  function  between  the  two  PSDs  as  the  difference  in  log  units  between  the 
candidate  and  reference  station  PSDs.  We  convert  the  data  from  decibels  into  an  amplitude  factor  that  converts  the 
candidate  data,  with  units  unknown,  to  the  units  of  the  reference  data  (e.g.,  acceleration).  This  amplitude  factor  is 
then  output  as  a  PAR  fde  (two  columns:  frequency  and  amplitude)  that  can  be  used  in  a  variety  of  seismic  software 
(e.g.,  SAC,  GSAC,  Antelope,  etc.)  to  remove  the  amplitude  response  for  the  reference  station  data. 
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Figure  1.  NOISETRAN  analysis  for  a  station  in  a  temporary  deployment  in  Israel.  (Left)  Candidate  (black) 
and  reference  (red)  station  noise  PSDs.  (Right)  The  transfer  function  (as  an  amplitude  factor)  to 
correct  the  candidate  PSD  to  the  reference  PSD. 


Figure  2  compares  the  PAR  obtained  from  NOISETRAN  (using  broadband  reference  station  EIL)  with  the  actual 
response  for  the  short-period  Sercel  L4  instrument.  NOISETRAN  has  successfully  recovered  the  true  instrument 
response  for  station  NS3  between  0. 1  and  3  Hz.  The  result  is  quite  remarkable,  considering  that  the  NS3  was  a 
temporary  deployment  buried  just  below  the  surface  while  EIL  is  a  pennanent  station  of  the  Israeli  Seismic 
Network.  The  differences  between  the  actual  response  and  PAR  at  high  frequencies  is  manageable,  better  than 
ROM,  and  could  be  tested  in  our  toolbox’s  next  stage,  which  includes  examining  signals  from  seismic  events  to 
further  refine  the  response.  These  results  suggest  that  NOISETRAN  will  work  for  short-period  stations  when 
corrected  using  a  nearby  broadband  station  with  similar  noise  conditions. 

Applying  NOISETRAN  to  USArray  Data 

In  the  second  test  of  the  Phase  II  project,  we  used  NOISETRAN  to  recover  the  response  of  both  transportable  and 
backbone  stations  of  the  USArray.  We  selected  a  one-hour  sample  of  nighttime  noise  from  1  February  and 
downloaded  the  data  from  IRIS  for  the  years  2009,  2010,  and  2011  to  use  as  input  for  NOISETRAN.  We  applied 
NOISETRAN  to  each  station  available  on  1  February  2009,  2010,  and  2011;  obtained  PARs  for  each;  and  analyzed 
errors  relative  to  the  known  responses.  Figure  3  shows  the  station  locations  for  the  years  201 1  (left)  and  2009  (right). 
Results  for  the  2011  dataset,  including  error  analysis,  are  presented  within.  As  an  example  of  the  NOISETRAN 
procedure.  Figure  4  shows  PSD  plots  for  a  candidate  station  (N38A)  and  a  reference  station  (Q34A).  These  stations 
are  separated  by  440  km.  The  difference  between  these  two  noise  spectra  in  dB  was  converted  into  an  amplitude 
factor.  Then,  a  PAR  was  generated  and  compared  to  the  known  frequency  response  of  station  N38A  (Figure  5).  This 
is  one  of  our  best  examples  of  the  agreement  between  the  NOISETRAN  estimate  and  the  true  response  for  a  station. 
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Cotnpanng  Pseudo  «xl  True  Resportses 


Fignre  2.  PAR  for  NS3  (±2x  Error)  from 

NOISETRAN  and  the  trne  velocity 
response  (hlack  line). 


In  our  tests  of  the  data  from  201 1,  we  performed  four  separate  analyses,  in  each  case  with  a  different  reference 
station,  in  order  to  test  the  range  of  performance  with  different  reference-candidate  station  geometries.  The  reference 
stations  selected  were  Q34A  (central),  BGNE  (off-center),  BLA  (eastern  US),  and  R1  lA  (western  US).  Station  BLA 
is  a  US  National  Seismographic  Network  (USNSN)  station,  while  the  others  are  from  the  USArray.  We  compared 
the  resulting  PAR  fdes  obtained  from  NOISETRAN  with  the  actual  station  responses  for  the  stations  shown  in 
Figure  5.  To  evaluate  the  method’s  performance,  we  calculated  mean  response  error  between  the  actual  response 
(^true)  and  the  PAR  (^par)  across  a  range  of  A  frequencies  (oji,  CO2,  . .  .con)  as 


er=^I]J:;i(loglO(dt^„e(a)))  -  loglO(dpAR(w)))-  (2) 

For  frequencies  between  0.01  and  10  Hz,  the  most  centrally  located  station  (Q34A)  provided  a  significantly  smaller 
amplitude  response  error  (Table  1)  for  the  mean  of  USArray  station  errors,  while  station  BLA  (USNSN)  provided  a 
significantly  smaller  error  for  the  mean  of  USNSN  station  errors.  Additional  error  analysis  for  reference  station 
Q34A  was  performed  to  evaluate  the  accuracy  of  PARs  at  using  lower- frequency  bands,  as  shown  in  Figure  6  and 
summarized  in  Table  1. 


USAnay  and  US  NationaJ  Slahona,  1  Feb  201 1 


USAfray  and  US  Nahorui  SlaMtm.  1  Feb  2006 


Figure  3.  USArray  stations  (red)  and  USNSN  stations  (blue)  on  1  February  2011  (left),  and  1  February  2009 
(right). 
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Noise  Plots 


Acceleration  Pseudo-Response 
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Fignre  4.  NOISETRAN  analysis  on  USArray  data.  (Left)  Candidate  N38A  (hlack)  and  reference  (red)  station 
Q34A  noise  PSDs.  (Right)  The  transfer  fnnction  (as  an  amplitnde  factor)  to  correct  the  candidate 
PSD  to  the  reference  PSD. 


Comparing  Pseudo  and  True  Responses 


Fignre  5.  PAR  for  N38A  (±2xError)  from  NOISETRAN  and  the  trne  velocity  response  (hlack  line) 


Table  1.  PAR  error  analysis  (2011) 


Reference 

Station 

Freqnency 
Band  (Hz) 

Mean  Error 
USArray 

Standard 

Deviation 

Mean  Error 
USNSN 

Standard  Deviation 

BONE 

0.01  -  10 

0.7390 

0.4612 

1.1400 

0.6557 

RllA 

0.01  -  10 

-1.3713 

0.4615 

-0.9609 

0.6606 

BLA 

0.01  -  10 

-0.5023 

0.4542 

-0.0968 

0.6622 

Q34A 

0.01  -  10 

-0.0554 

0.4618 

0.3476 

0.6557 

Q34A 

0.1  -  1.0 

-0.0857 

0.2282 

-0.0363 

0.5045 

Q34A 

0.01-0.1 

-0.0506 

0.1438 

-0.0348 

0.4701 
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We  plotted  the  PAR  errors  for  eaeh  station  versus  distanee  from  referenee  station  Q34A  for  the  USArray  stations 
(Figure  6)  and  the  USNSN  stations.  For  stations  within  200  km  of  the  referenee  station,  the  response  error  for  all 
frequeneies  is  typieally  within  less  than  ±0.3  log  units  (e.g.,  faetor  of  ±2  differenee  between  the  aetual  and  pseudo 
responses).  A  low-frequeney  band  (0.0 1-0. 1  Flz)  gives  even  tighter  error  bounds  that  extend  to  distanees  of  800  km 
or  more  from  the  referenee  station.  The  seatter  in  the  error  estimates  for  all  frequeneies  inereases  to  ±1.0  log  units 
(e.g.,  order  of  magnitude)  at  distanees  >  200  km;  however,  mueh  of  that  error  is  at  frequeneies  at  whieh  eultural 
noise  differenees  between  the  referenee  and  eandidate  stations  are  magnified. 


Mean  Error,  USArray  Stations  (0.01-10  Hz)  Mean  Error,  National  Network  (0.01-10  Hz) 
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Figure  6.  Mean  response  error  (red  crosses)  and  standard  deviation  (error  bars)  in  different  frequency  bands 
at  each  station,  plotted  versus  distance  from  reference  station  Q34A.  Top-right  shows  results  for 
USNSN  stations.  Top-left  and  bottom  plots  show  errors  for  USArray  stations  with  differing 
frequency  bands,  as  indicated. 
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Next,  we  plotted  the  mean  error  at  each  station  for  select  frequencies  in  map  view  (Figure  7)  with  an  overlay  of  US 
Geological  Survey  (USGS)  Generalized  Geologic  Map  of  the  United  States  (http://pubs.usgs.gov/atlas/geologic/). 
The  USGS  map  depicts  the  geology  of  the  bedrock  that  lies  at  or  near  the  land  surface,  but  not  the  distribution  of 
surficial  materials  such  as  soils,  alluvium,  and  glacial  deposits.  The  units  are  defined  by  sedimentary,  volcanic, 
plutonic,  or  metamorphic  rock  types  and  by  their  geologic  age.  For  this  exercise,  we  plot  only  the  unit  boundaries  as 
a  guide.  It  appears  that  there  are  correlations  between  the  amplitude  response  error  and  geology.  We  plan  to  analyze 
these  plots  further  to  better  understand  the  contribution  of  geologic  variations  to  response  error.  We  also  plan  to 
investigate  other  possible  error  sources  (e.g.,  topography  and  population  centers)  that  may  exist.  Any  correlations 
identified  can  help  us  to  make  more -appropriate  reference  station  selections  when  there  is  more  than  one  station 
under  consideration.  Additional  data  will  be  examined  to  determine  whether  observed  correlations  are  consistent  and 
can  be  used  to  develop  a  suitable  error  model  or,  alternatively,  a  correction  to  the  pseudo-amplitude  responses. 


True  vs  Pseudo  Amplitude  at  0.01  Hz 
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True  vs  Pseudo  Amplitude  at  0.05  Hz 


True  vs  Pseudo  Amplitude  at  0.1  Hz 
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Figure  7.  Mean  response  error  at  each  station  for  data  on  1  Feb  2011.  Color  scale  is  the  same  across  all  plots. 
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Testing  NOISETRAN  PARs:  The  11  March  Tohokn  Earthqnake 

To  further  test  the  PAR  files  generated  by  NOISETRAN,  we  applied  them  to  the  USArray  dataset  for  the  M  =  9.0 
1 1  March  201 1  Tohoku  earthquake.  PAR  files  were  created  for  each  station  in  the  March  2011  TA  network  (-430 
stations).  Station  Q34A  was  chosen  as  the  reference. 

In  order  to  define  a  consistent  P-arrival  for  the  complex  recordings  of  the  Tohoku  earthquake,  the  program  TauP 
was  used  to  calculate  the  P-wave  arrival  time  using  the  velocity  model  AK135  (Kennett  et  al.,  1995).  Waveforms 
were  windowed  approximately  600-900  sec  after  the  origin  time  of  the  earthquake  (3/1 1/2011  05:45:23  UTC)  to 
account  for  moveout  across  the  array.  The  PAR  files  were  converted  to  pseudo  phase-amplitude  response  files  for 
each  station  and  used  in  seismic  analysis  code  (SAC)  to  deconvolve  the  signal  in  the  time-domain.  The  same 
procedure  was  followed  for  the  true  response  files  (e.g.,  RESP  files  from  the  IRIS  SEED  volumes),  resulting  in 
“true”  and  PAR  corrected  waveforms. 

We  slightly  redefine  the  estimated  response  error,  as 

log lo(Ap  true)  “  loglo(Ap  p^;?),  (3) 

where  true  and  par  are  the  log  of  the  maximum  P-wave  amplitudes  in  the  analysis  window  for  the  true  and 

NOISETRAN-estimated  responses,  respectively.  The  response  error  is  small  for  the  unfiltered  data,  with  an  average 
-0.2  unit  bias;  however,  there  are  a  few  stations  that  are  clipping,  which  could  be  causing  this  low  error  estimate.  For 
bandpass  filtered  data  between  0.08  and  0.5  Hz,  the  error  increases  slightly  to  -0.3  log  units  but  increases  to  an 
average  -1.0  log  units  (e.g.,  order  of  magnitude)  with  a  0. 5-2.0  Hz  band  pass  filter.  The  bias  is  consistently  negative, 
meaning  that  the  PAR  files  are  producing  slightly  larger  P-wave  amplitudes  than  the  true  response.  The 
geographical  distribution  of  the  response  error  across  the  TA  network  using  reference  station  Q34A  is  shown  in  the 
maps  in  Figure  8  for  the  0.08-0.5  Hz  and  0. 5-2.0  Hz  filter  bands.  There  are  clear  geographic  patterns  observed  in 
the  response  error  that  need  to  be  further  studied  to  understand  the  performance  of  our  method.  We  are  surprised  by 
the  large  response  error  for  periods  near  1  sec,  and  this  might  suggest  poor  performance  of  the  NOISETRAN 
method  at  higher  frequencies  or  poor  choice  of  a  reference  station.  We  will  continue  to  examine  this  dataset  in  the 
coming  months  of  the  project,  including  using  synthetic  analysis  of  the  event  to  iterate  to  an  improved  response 
function. 
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Fignre  8.  Geographic  distrihntion  of  response  error  for  Tohnkn  /’-waves  filtered  at  0.08-0.5  Hz  (top)  and 
0. 5-2.0  Hz  (hottom)  nsing  reference  station  Q34A. 


CONCLUSIONS  AND  RECOMMENDATIONS 


After  the  first  nine  months  of  the  Phase  II  projeet,  we  have  made  signifieant  progress  in  the  applieation  of  SACPSD 
and  the  development  of  NOISETRAN,  whieh  are  the  major  eomponents  of  our  response  reeovery  toolbox.  The 
progress  ineludes  the  following: 

•  We  now  understand  the  distanee  requirements  for  the  referenee  stations  (e.g.  ,<  200  km)  and  the  magnitude 
of  the  response  error  that  ean  be  expeeted  aeross  a  frequeney  band  of  0.01-10  Hz.  We  also  now  have  a 
better  understanding  of  geologieal  and  other  effeets  on  the  noise  eomparisons  between  the  referenee  and 
eandidate  stations.  For  periods  below  1  Hz,  the  referenee  station  seleetion  ean  be  extended  beyond  200  km. 

•  We  have  also  developed  a  method  for  eonverting  the  PAR  into  an  amplitude  file  that  ean  be  used  in  SAC, 
GSAC,  Antelope,  or  other  programs  for  easy  seismometer  response  deeonvolution.  This  allows  us  to 
eorreet  real  signals  (e.g.,  P-waves  from  the  Tohuku  earthquake)  using  the  PAR  flinetions  estimated  from 
noise  analysis. 
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•  We  have  tested  new  ways  for  metadata  evaluation  and  response  error  deteetion. 

•  We  are  eurrently  working  again  on  a  graphieal  user  interfaee  for  the  toolbox. 

During  the  final  year  of  the  projeet,  we  will  paekage  NOISETRAN  and  additional  seismie  tools  into  an  easy-to-use 
graphieal  interfaee,  eonduet  alpha  testing  of  the  toolbox,  and  then  eonduet  beta  testing  at  Los  Alamos  National 
Laboratory.  We  also  need  to  examine  possible  phase  reeovery  of  the  seismometer.  Based  on  our  initial  presentation 
of  these  results  to  other  seismie  researehers,  we  believe  this  toolbox  will  have  an  important  impaet  in  many 
seismologieal  eommunities  where  the  eharaeteristies  of  the  reeording  system  are  either  unknown  or  possibly 
erroneous. 
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